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Abstract 

Near conformal dynamics is employed in different extensions of the standard model of particle inter- 
actions as well as in cosmology. Many of its interesting properties are either conjectured or determined 
using model computations. We introduce a relevant four dimensional gauge theory template allowing us 
to investigate such dynamics perturbatively. The gauge theory we consider is quantum chromodynamics 
with the addition of a meson-like scalar degree of freedom as well as an adjoint Weyl fermion. 

At the two-loop level, and in the Veneziano limit, we firmly establish the existence of several fixed points 
of which one is all directions stable in the infrared. An interesting feature of the model is that this fixed 
point is lost, within the perturbatively trustable regime, by merging with another fixed point when varying 
the number of quark flavors. We show the emergence of the Miransky scaling and determine its properties. 
We are also able to determine the walking region of the theory which turns out to be, at large number of 
colors, about 12% of the conformal window. Furthermore, we determine highly relevant quantities for near 
conformal dynamics such as the anomalous dimension of the fermion masses. 
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I. INTRODUCTION 



Determining the phase structure of gauge theories of fundamental interactions is crucial in 
order to be able to select relevant extensions of the standard model of particle interactions for 
particle physics, dark matter [1], cosmology [2, 3] and quantum gravity [4]. 

In particular it is interesting to elucidate the perturbative and non-perturbative physics of 
non-Abelian gauge theories with and without fundamental scalars as function of the number of 
flavors, colors and matter representation. 

The relevant points to address are: 

i) Determine the regions of conformality of a given gauge theory as function of the couplings, 
number of flavors and matter representation. 

ii) Determine whether near the conformal boundary between a conformally broken and con- 
formally restored symmetry the theory walks or jumps? Walking dynamics [5-8] requires 
that the physical observables of the theory, in the broken phase, feel the presence of the 
conformal phase. This must happen at least in a region of the parameter space of the theory 
near the phase transition. Jumping dynamics [9] refers instead to the logical possibility that 
the theory in the broken phase displays no signs of a nearby, in the parameter space, con- 
formal region. Jumping dynamics is associated to a first order conformal phase transition 
while walking dynamics to a second order or higher transition. 

iii) Jumping and walking dynamics could be realized physically in different ways. For example 
walking arises naturally when an ultraviolet fixed point (UVFP) and an infrared one (IRFP) 
smoothly merge. The full beta function of the theory near this merger can be modeled by the 
product of two zeros. This mechanism also underlies the conjectured Miransky scaling [10- 
12] describing how the intrinsic scale of the theory vanishes as we approach the conformal 
boundary. On the other hand the ratio of two zeros in the beta function leads naturally to 
jumping dynamics [9]. The problem is therefore to establish which mechanism underlies 
either walking or jumping dynamics in a given gauge theory. 

iv) Once the kind of the conformal transition has been established one needs to determine the 
physical properties of the theory. These are, for example, the anomalous dimensions of the 
fermion mass operator and the particle spectrum. 



2 



v) If the conformal transition is of walking type another question arises: What is the actual 
range of the parameter space of the theory supporting walking, e.g. the walking window 
as function of the number of flavors and colors? 

To help answer some of the points above it is, therefore, relevant to construct calculable 
examples such as the two dimensional toy model investigated in [13]. Here, however, we are 
interested in investigating four dimensional gauge theories allowing us to address all, or part, of 
the questions above [14, 15] in perturbation theory 1 . Since we are not concerned with the hierarchy 
problem, this is obviously true also for the toy model studied in [13], we consider theories featuring 
also fundamental scalars besides fermionic matter as done in [14-16]. 

We start by introducing the gauge theory template [15] allowing us to investigate conformal 
dynamics and its breaking within perturbation theory. The gauge theory we consider is quantum 
chromodynamics (QCD) with the addition of a meson-like scalar degree of freedom as well as 
an adjoint Weyl fermion. This model has a large number of global symmetries in common with 
ordinary QCD and allows us, at the two-loop level, and in the Veneziano limit, to establish the 
existence of several fixed points of which one is all directions stable in the infrared. 

A relevant feature of the model is that the infrared stable fixed point can be lost, within the 
perturbatively trustable regime, by either the merging with a non-trivial ultraviolet fixed point 
when varying the underlying number of quark species or via spontaneous condensation 2 of a 
scalar field due to perturbative quantum corrections [17]. Interestingly these two mechanisms 
associated to the loss of an infrared fixed point can be investigated in great detail within the same 
gauge theory. The merger phenomenon [9, 18-22], as we shall see, leads to the expected slow 
running of the coupling constant towards the infrared also known as walking dynamics pioneered 
by Holdom [5, 6] and also investigated in [7, 8, 23, 24]. The fixed points merge when the number 
of flavors reaches a critical value. Due to the calculable nature of the model we can provide a clear 
realization of the Miransky scaling [10-12] and walking dynamics. 

The phenomenon of merging fixed points unveils itself here when going to the second order in 
perturbation theory. In fact, to the first order conformality can be lost, in our model, only because 
of the quantum corrections inducing spontaneous condensation of the scalar degrees of freedom 
[15]. Other quantities for walking dynamics such as the anomalous dimension of the fermion 
masses will also be determined precisely. 

1 Actually in [15] we showed that one can also address some non-perturbative aspects. 

2 We use the short-hand notation "condensation of the scalar field" to mean the spontaneous generation of its vacuum 
expectation value. 
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In section II we review the model and the one-loop results. In section III we setup the two-loop 
beta function analysis and discuss the fixed point structure. We also show that the IR stable fixed 
point merges with an UV fixed point as we decrease the number of flavors within the perturbative 
regime. We then move to section IV where we investigate in detail the physics of the merger. Here 
we determine precisely the conformal window as function of the number of flavors. Remarkably 
we are even able to determine the precise region of walking. We show that the Miransky scaling 
emerges naturally in this model. We finally compute that anomalous dimension of the mass 
operator both for the UV and IR non-trivial fixed points. We briefly conclude in section V while in 
the appendix A we provide the full beta function and anomalous dimension of the fermion mass 
operator of the theory at two loops and for any number of colors. 

II. QCD WITH ONE ADJOINT FERMION AND ELEMENTARY MESONS 

We want to investigate the dynamics of the following gauge- Yukawa theory: 
1 



£ = Tr 



-^F f,v F fI v + iA0A + Q0Q + d^dm + y H QHQ 
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- U!(Tr[H + H]) 2 - u 2 Tr(H + H) 2 . (1) 



The theory is, at the tree-level, conformally invariant. With "Tr" we indicate the trace operator over 
both color and flavor indices and is the covariant derivative respecting the gauge symmetries 
of the field on which it acts. The field content and symmetries are given in Table I, where the Dirac 
fermion is decomposed in its Weyl components Q = (q, q*). The gauge coupling constant will be 
denoted by g. This theory is QCD with one adjoint Weyl fermion and elementary mesons. It has 
been chosen since, as we shall see, it allows for a very interesting perturbative gauge dynamics, 
and for its close resemblance to QCD. 

TABLE I: Field content. The first three fields are Weyl spinors in the (k,0) representation of the Lorentz 
group. H is a complex scalar and G f , are the gauge fields. U(1)af is the extra anomaly free symmetry 
arising due to the presence of A. 



Fields 


[SU(N C )] 


SU(N f ) L 


SU(N f ) 


A 


Adj 


1 


1 


<7 


□ 


□ 


1 


? 


□ 


1 


□ 


H 


1 


□ 


□ 




Adj 


1 


1 



N c N f 
Nf-Nc _N £ 
N c N f 

N f 
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Throughout this paper we will work with the rescaled couplings which allow for a finite 
Veneziano limit of the theory, i.e. N c ,Nf — > oo while keeping x — Nf/N c fixed. The appropriate 
rescaled couplings read 

_ g 2 N c _ y 2 H N c _ UiNj _ U2 N f 

H ~ (471)2 ' UH ~ (471)2 ' 21 - (471)2 ' 22 " (471)2 ' 

If not otherwise stated, we will be working in the Veneziano limit. 

This model was introduced in [15], where a one loop analysis was performed and showed to 
yield several interesting results. In this paper, we start the two-loop analysis of the theory where, 
as we will demonstrate, new important features emerge. It is relevant for the completeness of this 
paper to summarize the most important features found at one loop, while referring the reader to 
[15] for the details of the analysis (see also [14] for another model featuring similar properties). 

i. The theory has a perturbative, infrared stable fixed point (IRFP) in all couplings near the 
asymptotically free boundary given by: 

x = 9/2 . (2) 

The perturbative nature of the fixed point (FP) is illustrated by showing the value of the 
couplings at the IRFP as a Taylor expansion around the asymptotically free boundary x to 
leading order: 

22 _ 4 _ -2 V19 + J2(8 + 3V19) _ _j + yi9 _ 
4 = ~ a H = ~ x )> z i = 27 ( x - z *2 = 27 ^ ~ ^ ' 

We note that the existence of a perturbative stable IRFP for a generic non-supersymmetric 
gauge theory featuring scalars is not automatic [25]. For instance, if we were to eliminate 
the adjoint Weyl fermion A and consider simply QCD with elementary mesons, such a 
perturbative fixed point would disappear. In Fig.l we show the dependence of the IRFP on 
x to better appreciate its perturbative nature. 

ii. Although there is a perturbatively stable IRFP not all the Renormalization Group (RG) 
trajectories will be able to end on it because of the presence of elementary mesons. In fact, 
the scalar potential can be unstable against a spontaneous chiral symmetry breaking phase 
transition induced by the Coleman- Weinberg (CW) phenomenon [17]. This phenomenon is 
not encoded in the perturbative running of the couplings. A critical boundary (separatrix) 
was discovered in the parameter space of the theory separating the stable region from the 
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unstable one. This is shown as a two-dimensional section of the parameter space in Fig. 2, 
where it is assumed that the gauge and Yukawa couplings have quasi-reached the FP. 

iii. In the stable region, i.e. where the CW does not lead to spontaneous symmetry breaking, 
the FP is reached and the theory exhibits large distance conformality. In the unstable 
region the effective potential develops a non-trivial minimum leading to chiral symmetry 
breaking to its maximal diagonal subgroup. The resulting perturbative spectrum consists 
of heavy fermions, heavy scalars, massless Nambu-Goldstone (NG) bosons and a pseudo- 
NG boson associated with spontaneous breaking of the dilatation symmetry. The gauge 
sector and A-field are not affected by the symmetry breaking and emerge as a low-energy 
super Yang-Mills theory of which the condensate is computable nonperturbatively. The 
scale of the condensate/confinement of this remaining low energy theory can be shown to 
be exponentially lighter than any other generated mass scale of the theory. 

iv. It has also been shown that the pseudo-NG boson, associated to dilatonic invariance, is the 
dilaton state of the theory at the quantum level. Its mass, to one loop order, was computed 
to be: 



Zl = 



r , ( l Azj\n(Az 2 ) - xa 2 \n(xa H y 

V^-x^Jexp 2 



4z 2 - xajj 



(3) 



where /Jo is the scale at which the (bare) Lagrangian is defined and it is approximately the 
mass scale of the other heavy states. It was shown that m 2 can be made much lighter than 
the heavy mass scale by tuning 4z 2 — xa^ «: an, z 2 , since the heavy states acquire tree level 
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FIG. 1: Fixed point solutions at one loop as a function of x = Nf/N c . 
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FIG. 2: Showing the separation of phases of the RG flow, as discussed in the text, for x = 4.4. The other 
couplings, dg and an, are imagined to have quasi-reached the fixed point for 

masses proportional to an and Z2- This tuning corresponds to RG trajectories closer and 
closer to the separatrix. A large scale separation between the different states of the theory is 
therefore obtained in the chirally broken region by tuning the theory closer to the separatrix. 

We have thus constructed an explicit perturbative realization of conformal and near conformal 
dynamics. However, this realization differs from the one in QCD, since it is perturbative and 
furthermore the tuning needed to go in and out of the conformal window is in the couplings 
instead of the number of quark flavors. Furthermore, the primary condensation phenomenon 
involves a scalar field rather than the quark condensate. As we shall see in the following sections, 
this situation becomes more interesting when two loop corrections are included. 

III. THE TWO LOOPS BETA FUNCTIONS AND EMERGENCE OF NEW FIXED POINTS 

It is interesting to go beyond one-loop since we know that higher order analyses often reveal 
several new important features [26-33] not present at the lowest order. In fact, by carrying 
the analysis one step further in the perturbative expansion we will discover that the IRFP can 
also disappear, within the perturbative regime, by merging with a new ultraviolet fixed point 
as we decrease the number of quark flavors. The mechanism behind the loss of conformality is 
very different from the case in which this occurs via scalar condensation induced by quantum 
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corrections. In this case, it is reasonable to expect the quarks to condense first and induce the 
scalar condensate. Of course, we have to work in the region of the parameter space where the CW 
potential does not lead to a scalar condensate. We have checked, by an explicit computation, that 
such a region of couplings exists. 

We now therefore present the beta functions and their fixed point (FP) structure at two loops. 
The two-loop beta functions have been derived using the formalism of [34-37]. In particular, the 
parametrization in [38] was used to handle the scalar sector. The beta functions read: 



p Bg = - - a 2 [ll - 21 - 2x + (34 - 16£ - 13x) a g + 3x 2 a H ] (4) 

\ , , 20(x + f)-203 7 x{x + \2) 7 J , x 
pa H = an 2(x + - 6a g + (8x + 5)a g aH -\ a g - Sxz 2 an ^ <% + ^2 (^) 

p Zl = 4(z* + 4zjZ 2 + 3z^ + z x a H ) 

(6) 

+ 2 {5zia g a H + 2x 2 a\ i + x(4z 2 - 3z 1 )a 2 j - Az\a H - \2z\a H - l6ziZ 2 a H - 4&z\ - 20z 1 z 2 ^ 

fi Z2 = 2 (2z 2 ah + 4z^ - xa 2 ^ + 2 (5z 2 a g a H - 2xa g a 2 H + 2x 2 a 3 H - 3xz 2 a 2 H - Sz\a H - 12z^) (7) 

The parameter i is the number of A fields of the theory. In our case i = 1. We have kept t as an 
external parameter, such that the results can be used for future research e.g. the model of QCD and 
mesons (£ - 0). Notice that the double trace coupling z\ decouples from the beta functions of the 
remaining couplings a g , an, z 2 in the Veneziano limit. Moreover, j$(z\) is a quadratic equation in Z\. 
Both of these properties actually hold to all orders in perturbation theory in the Veneziano limit, as 
was shown in [39]. In Appendix A we provide the full Nj dependent beta functions to two loops. 
We have determined these expressions by implementing the above mentioned formalism in two 
independent ways for cross checking and have furthermore verified them by comparing with the 
partial two loop results obtained in [40] for a similar model without gauge interactions and the A 
field. We have also checked that the expressions are consistent with the Callan-Symanzik equation 
for the two loop effective potential [41]. 

The system of beta functions has a rich fixed point structure which we will briefly discuss. 
Because of the decoupling and polynomial structure of |S(zi) mentioned above, we know that it 
will always either have two real or no (i.e. two complex) FP solutions. For the remaining system 
of beta functions, we have found, depending on the value of x, at most six real solutions. One of 
these corresponds to the perturbative IRFP discussed in Sec.II. We focus our attention on the two 
loop fate of this IRFP as function of the number of flavors, i.e. x. To help the reader we provide 
its explicit two loop dependence on the number of flavors in Fig.3. It is clear from the figure that 
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FIG. 3: Merging fixed points at two loops. 

the IRFP merges at a critical number of flavors with a new ultraviolet (UV)FP appearing at two 
loops. Interestingly the UVFP appears to be perturbative in the same range of X in which the IRFP 
exists. Moreover, the merger occurs at x* = 4.42028 which is quite close to the asymptotically free 
boundary for which the perturbative expansion is valid. It is therefore natural to expect that higher 
loop effects will not qualitatively affect this picture. This means that we discovered a conformal 
window as function of number of flavors accessible in perturbation theory. 

IV. PERTURBATIVE MIRANSKY SCALING, WALKING DYNAMICS AND CRITICAL EXPONENTS 

Having investigated the fixed point structure of the theory at two loops we can turn our 
attention to the primary questions we address here. What kind of near-conformal dynamics 
appears as we change the number of flavors? Does the theory jump or walk? We will demonstrate 
that: a) The investigated gauge theory displays a smooth phase transition as a function of the 
number of flavors when crossing the lower boundary of the conformal window; b) There is a well 
defined region of walking dynamics; c) We can reliably compute the relevant critical exponents of 
the theory such as the anomalous dimension of the fermion mass operator in the entire conformal 
window. 

We have already discussed that the IRFP as function of the number of flavors merges with the 
UV one, and argued that this happened within the perturbative regime of the theory. It is therefore 
timely to start a careful perturbative investigation of the physical properties of the theory around 
the lower boundary of the conformal window delimited by the merging of the IRFP and UVFP 
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a 

FIG. 4: Beta functions for a = a gr an,zi,Z2 and for x = 4.425. Each beta function features two 
nontrivial FPs, one in the UV, which is unstable under small changes in the values of the four couplings, 
and one stable in the IR. 

occurring at x* = 4.42028 corresponding to a critical number of flavors N* = x*N c . Asymptotic 
freedom is lost for Nf - xN c - 9N c /2. We will also consider what happens when we lower x past 
the point of the merger in order to demonstrate the emergence of walking dynamics. 

Within the conformal window, i.e. x* < x < x the physics is controlled by the nontrivial IRFP 
point when starting the flow in the ultraviolet near the trivial UVFP. In Fig. 4 we show the different 
fl a (a) for a - a g ,ciH,zi,Z2 at x - 4.425 obtained by solving the RG equations numerically 3 . 

When x is only slightly lower than x* , we find that this situation is mostly unchanged, since, 
as we shall demonstrate, the dynamics is controlled by the quasi FPs. The beta function jS a (a g ) is 
shown in Fig. 5 for different values of X near the merger. 

Below the merger the theory must develop a physical scale A and we introduce therefore the 
renormalization group time t = log ^ it takes for the theory to flow from a UV scale fio to the infrared 
scale A. 

It is clear from Fig. 5 that a local maximum of is present at a value of the coupling close to 
the merger FP value. The range of values of x < x* for which jS a features a local maximum away 
from the trivial FP defines concretely the walking dynamics regime. Below the value of x where 
the local maximum disappears the theory no longer walks but runs. From Fig. 5 it is evident that 

3 Specifically each of the four betas is composed of three parametric plots of § a {^i),a{[i), determined by choosing 
different initial conditions, respectively one close to the origin, and two close to the UV FP but on different branches 
of p a , and for a large range of the renormalization scale /.(. 
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FIG. 5: The beta function jS flg (fl ? ) plotted for x = 4.42,4.4167,4.4133,4.41. The local maximum of jS % 
disappears for x 4.41. 



this happens for x « 4.41. We have thus managed to estimate the walking window of the theory: 

< x < 4.41 Running Dynamics 

4.41 < x < 4.42 Walking Dynamics 

4.42 < x < 4.5 Conformal Window 
4.5 < x IR free 

The phase diagram including the conformal window (light blue shaded region) and walking 
region (pink shaded region) as a function of the number of flavors and colors is reported in Fig. 6, 
taking into account the full dependence on the number of colors N c . Deviations from the large N c 
limit (dashed line) are sizable at small N c as it can be seen from the figure. 
We can then determine A in the walking regime as follows 

^A-"., ra g (A) da Q 



log 



A 



(8) 



Mo J fIO f J ag(flo) Pa g (a 8 ,h,h,h) 

where ki with i — 1, 2, 3 are the integration constants for the other couplings solved as a function 
of flg. Within the walking region the theory will spend most of the RG time near the quasi-FP a* , 
corresponding to a local maximum of p ag . It is therefore natural to pick fj.o such that flg(/io) = a *g 
[42]. In the deep infrared the a s coupling diverges defining A. We have computed the last integral 
in Eq. 8 numerically, which is shown in Fig. 7. It is interesting to compare the numerical result 
with the leading order analytic contribution of |S fl? to the integral in Eq. 8, which is determined by 



11 



expanding jS a around its local maximum: 




with 

The quantities jSo,2 depend on x and, of course, also on the Jq integration constants. We evaluated 
them for 4.41201 < x < 4.42021 in increments in x of 10~ 4 and define fio,2(x) as the best fit fourth 
order polynomials in x—x* to the respective data points 4 . The resulting leading order contribution 
in x we obtain for the RG time is 

log — « , a = 2.99 x 10" 2 . (11) 

Mo Ci V** - x 

The result above shows that the model realizes the Miransky scaling [11]. The Miransky scaling 
function in Eq. 11 is plotted in Fig. 7 together with the precise numerical evaluation of the integral 
in Eq. 8 for the interval 4.41201 < x < 4.42021, which covers most of the walking region in x. It is 
clear from Fig. 7 that the Miransky scaling approximation becomes increasingly accurate as x gets 
closer to x*. This is an expected result, since for x — > (x*)~ higher order corrections to the inverse 
parabolic approximation of fi~ 1 become negligible. 

4.50 
4.45 

N f 



4.35 

4 6 8 10 12 14 16 18 

N c 

FIG. 6: Phase diagram including the conformal window (light blue shaded region) and walking region 
(pink shaded region) of the theory as function of the number of flavors and colors. The dashed line 
corresponds to the large N c limit of the lower boundary of the conformal window. 



4 Higher order polynomials do not change our results. 
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FIG. 7: The Miransky scaling function (red dashed) in Eq. 11 plotted together with the numerically 
precise evaluation of the integral in Eq. 8 (solid blue) for the interval 4.41201 < x < 4.42021. The 
Miransky scaling approximation becomes increasingly accurate as X gets closer to x* . 

So far we have shown that the model possesses a smooth conformal phase transition when 
crossing the lower boundary of the conformal window as function of the number of flavors. 
At the FPs anomalous dimensions are expected to be physical quantities identifiable with the 
critical exponents of the theory. These are important quantities governing the scaling properties 
of relevant correlators and therefore deserve to be investigated. A particularly interesting one is 
the anomalous dimension of the fermion mass operator. This quantity evaluated at the merger is 
very relevant for walking dynamics 

dlogm 

ym ^"di^g7' A # = 3-7«- (12) 
The direct contributions to y m at the one and two-loop order are [37]: 

y^' = 3a g - xuh , (13) 
y J = — (96xa g a H + (203 - 20x - 20£)a 2 g + 3x(x + 6)a 2 H ) . (14) 

As in the expressions for the beta functions t denotes the number of A fields. For completeness, 
the full two-loop expression for the anomalous dimension, not truncated to the Veneziano limit, is 
provided in appendix A. Notice that the quartic interactions do not contribute to the anomalous 
dimension at two-loop order. In Fig. 8 the red and black lines corresponds to y m = yjp + y^ at the 
IRFP (red) and UVFP (black) given in Fig. 3. The physical dimension of the fermion mass operator 
at the merger reads « 2.92. In the walking regime the physics is dominated by the nearby FPs 
and therefore the relevant anomalous dimension is y m erger w 0.08. 
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FIG. 8: Anomalous dimension of the fermion mass at the UVFP (upper black) and IRFP (lower red). 

Remarkably this model features all the salient properties of walking dynamics. Of course, 
larger values of the anomalous dimensions can arise for theories which are strongly coupled at 
the phase transition point. Nevertheless it is interesting to have shown that walking dynamics is 
actually realized in asymptotically free four-dimensional gauge theories. 

V. CONCLUSIONS 

We introduced a relevant four dimensional gauge theory template which has allowed us to 
determine its conformal and walking windows as functions of the number of flavors within 
perturbation theory. We have shown that conformality is lost due to the merging of an IR and 
UV fixed points leading to Miransky scaling. We were able to define the walking regime of the 
theory and compute the dependence on the number of flavors of the anomalous dimension of the 
fermion mass. 

Walking dynamics is required as an essential feature for extensions of the standard model 
featuring new asymptotically free gauge theories. However, many of its required properties are 
either conjectured or investigated using model computations. Here we have first introduced a 
relevant asymptotically free gauge theory featuring Nf flavors in the fundamental representation 
of an SU(N C ) gauge group, like QCD, and then added a Weyl fermion in the adjoint together with 
meson-like fundamental scalars. 

At the two-loop level, and in the Veneziano limit, we discovered the existence of several fixed 
points of which one is all directions stable in the infrared. We showed that this fixed point is lost, 
within the perturbatively trustable regime, by merging with another ultraviolet fixed point when 
varying the number of quark flavors. We naturally recover the Miransky scaling and determine 
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its properties. The theory is manifestly of walking type and the theory allows us to determine 
the associated walking window in the number of colors and flavors. We find that the walking 
window is about 12% of the conformal window at large N. Furthermore, we determined also other 
interesting quantities for walking dynamics such as the anomalous dimension of the fermion mass. 

Appendix A: The full two-loop ^-functions and anomalous dimension 

For completeness, we provide the jS-functions of the theory and the anomalous dimension of 
the fermion mass to the full two loop order. The parameter ( denotes the number of A fields. 
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